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Предложена математическая модель нагрузки импульсного наземного источника сейсмических волн, 
позволяющая выбирать параметры основных элементов его конструктивной схемы при воздействии на грунт 
с превышением предела его упругости. 
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Введение. Сейсморазведка является одним из основных методов поиска полезных ископаемых. 
Эффективность сейсморазведочных работ определяется в значительной мере характеристиками 
источников сейсмических сигналов. Очевидным способом генерации сейсмических волн с поверх- 
ности грунта является деформация ее источником мощных механических воздействий (импульс- 
ных или вибрационных). 

В конструкциях наземных импульсных сейсмоисточников, как правило, есть помещенная 
на поверхность грунта плита — излучатель, — над которой расположен массивный груз. Податли- 
вость излучателя и груза значительно меньше податливости грунта. Для импульсного воздействия 
через плиту на грунт груз сбрасывают на нее с заданной высоты (первый способ) или в сейсмо- 
источнике срабатывает специальный мощный импульсный привод (в настоящее время использу- 
ется в основном электромеханический), отталкивающий груз от плиты (второй способ) [1]. 

Большое значение при проектировании автономного импульсного сейсмоисточника 
с ограниченным запасом энергоносителя имеет эффективность преобразования энергии привода 
или падающего груза в энергию сейсмических волн. Известно, что лишь часть механической энер- 
гии, переданной излучателю, преобразуется в целевые сейсмические волны — объемные продоль- 
ные или поперечные. Часть ее преобразуется в поверхностные волны, являющиеся при сейсмо- 
разведке помехами. Распределение энергии генерируемых волн у современных сейсмоисточников 
с перпендикулярным поверхности воздействием на грунт такое: продольные - 7-9%, попереч- 
ные - 18-25%, поверхностные - 68-75% [2, 3]. Столь нерациональное использование механиче- 
ской энергии в немалой степени обусловлено отсутствием возможности оптимизировать парамет- 
ры воздействия сейсмоисточника на грунт на этапе проектирования, а подбор размеров излуча- 
теля изготовленного серийного сейсмоисточника экономически нецелесообразен. В сложившейся 
ситуации ведущие сейсморазведочные работы организации вынуждены использовать сейсмои- 
сточники, имеющие запас по мощности и генерирующие волны в широком диапазоне частот. 

Расчетные модели взаимодействия вибрационного источника с полупространством 
грунта, позволяющие оптимизировать параметры воздействия на грунт, описаны [4, 5], а для им- 
пульсных сейсмоисточников адекватных указанных моделей нет, что является проблемой. Моде- 
лирование и исследование процесса взаимодействия наземного импульсного сейсмоисточника с 
нагрузкой — грунтом — актуальны и имеют прикладную и научную значимость при постановке но- 
вых задач сейсморазведки и при проектировании сейсмоисточников. 
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Цель исследования — разработка простой и адекватной математической модели 

нагрузки импульсного наземного сейсмоисточника в соответствии с поставленными задачами 
сейсморазведки, пригодной для оценки соответствия параметров формируемого механического 
импульса (величины, формы и длительности) размерам излучателя. 
Методики моделирования процесса генерирования сейсмических волн механическим 
воздействием. Общее решение задачи о распространении в однородном упругом полупростран- 
стве произвольного воздействия, приложенного к свободной поверхности через излучатель огра- 
ниченной площади, известно как задача Лэмба. Согласно этому решению, при импульсном воз- 
действии на полупространство в нем возникают и распространяются волны четырех типов: объ- 
емные продольные и поперечные, поверхностные и обменные, обусловленные взаимными преоб- 
разованиями волн разного типа на границах областей их существования [4]. Эти волны распро- 
страняются от излучателя по различным секторам и с различными скоростями: самая большая 
скорость у продольных волн, самая маленькая - у поверхностных. 

Объемная продольная волна в грунте, генерируемая расположенным на нем излучателем, 
распространяется в секторе, ось которого совпадает с перпендикуляром к поверхности в центре 
излучателя. В пределах этого сектора для математического описания процесса генерации сейсми- 
ческой волны можно использовать модель «сферического источника» [4] (реализуемого взрывом 
в скважине либо пульсирующим в воде газовым пузырем). Эта модель дает математическую трак- 
товку нескольких простых физических принципов, что позволяет на их основе формулировать 
практические рекомендации. Для объемной продольной волны характерно следующее: 

— имеет на любом расстоянии от источника определяющийся характером поля смещений 
на границе источника так называемый «график колебаний». Характер поля смещений источника 
представляет собой сравнительно быстро затухающие (в течение трех-четырех периодов) коле- 
бания; 

— в пределах зоны не более одной длины волны от излучателя, кроме смещений грунта 
колебательного характера, обусловленных прохождением волн сжатия и разрежения, присутст- 
вуют также смещения вынужденного характера — апериодические. «График колебаний» в этой 
зоне содержит составляющие, пропорциональные и характеру поля смещений (с коэффициентом 


у "?, где г — расстояние от излучателя) и первой производной от него (с коэффициентом Ш) 


[4]. На значительном удалении от источника (более одной длины волны) из-за затухания вынуж- 
денной составляющей «график колебаний» приобретает практически неизменную форму, опре- 
деляющуюся в основном первой производной от характера поля смещений (скоростью смещения). 
Регистрация сигналов в сейсморазведке осуществляется датчиками скорости; 

— так как сила воздействия импульсного сейсмоисточника на грунт изменяется во време- 
ни, то в «графике колебаний» на значительном удалении от излучателя можно выделить две 
составляющие. Одна определяется размерами излучателя и волновыми свойствами грунта, а вто- 
рая — скоростью движения излучателя. 

Однако в модели «сферического источника», как и во многих других математических мо- 
делях и решениях задач теории упругих колебаний, используется однородная среда. С использо- 
ванием этой модели сложно установить соответствие характеристик излучаемого импульсным 
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сейсмоисточником сигнала параметрам механического импульса, а также размерам излучателя, 
так как не учитывается множество следующих важных факторов, например: 

— невысокие (менее 5`10? Па) значения предела упругости в поверхностном слое грунта — 
в так называемой зоне малых скоростей (3МС) - и уменьшение модулей упругости (Юнга, Ламе) 
при уменьшении скорости воздействия; 

— многокомпонентная неоднородная структура грунта (газ, жидкость, твердые фракции 
различных форм и размеров); 

— свойство грунта уплотняться в ходе каждого сжатия и разрыхляться при снятии с него 
давления, т.е. изменчивость его структуры -— «скелета» - ит. д. 

Таким образом, использование структурных распределенных математических моделей на- 
грузки импульсного сейсмоисточника, в которых объектами исследования являются поля физиче- 
ских величин, ограничено. Число совместно исследуемых различных сред (число деталей, мате- 
риалов, фаз агрегатного состояния) в практически используемых моделях такого уровня не может 
быть большим ввиду сложностей вычислительного характера. Кроме того, для исследования 
только лишь динамических характеристик сейсмоисточника такая модель избыточна. Для этой 
цели может оказаться достаточной модель грунта с сосредоточенными параметрами. Для ее по- 
строения необходимо провести ранжирование свойств грунта как нагрузки сейсмоисточника и 
включать в ее модель наиболее значимые свойства. 

Свойства поверхностного слоя грунта. Основным ограничивающим фактором процесса гене- 
рации сейсмических волн в грунте, который обязательно должен быть учтен в модели, является 
предел упругой деформации с, грунта. Зависящие от величины с, предельные значения пара- 


метров воздействия сейсмоисточника на грунт определяются по известным формулам: 
— максимальная скорость смещения частиц грунта в продольной волне упругих деформа- 


ций, задающая предел скорости У движения излучателя [6]: 


У, =6,/(ру,, (1) 


где ру, — волновое сопротивление грунта, Па’с/м; р — плотность грунта, кг/м; у, — скорость 


шах 


распространения продольных сейсмических волн в грунтовом полупространстве, м/с: 





Здесь Е - модуль Юнга, Па; и - коэффициент Пуассона. При скорости движения излучателя 


большей у „. грунт под ним подвергается пластическим деформациям; 


Хх 

— максимальное время действия деформаций сжатия определяется временем, после кото- 
рого избыточное давление в волне напряжений в грунте, контактирующем с движущимся вниз 
излучателем, исчезает [7]: 


:„„” Б/у, , (2) 
где р - диаметр излучателя (для излучающей поверхности круглой формы), м: 
р=2./5./п, 
5, — площадь излучающей поверхности, м": 
Вы О (3) 


Здесь Р’„. - максимальное значение развиваемого сейсмоисточником усилия Ё воздействия на 


грунт. При превышении Р„ также возникают пластические деформации грунта; 
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— минимальное время действия деформаций [6]: 
ан = Х пах У = ХпахРУ, /с, ' (4) 
где х„‚ — максимальное значение упругих деформаций грунта под излучателем, м. Если пренеб- 


речь неравномерностью нагрузки по площади излучающей поверхности, то получим: 
хх ®Ро,/Е, 


в =РУ,Б/Е. 
Известно, что предел упругой деформации для неуплотненных поверхностных слоев грун- 
та весьма мал: о, =(2-5)-10° Па, и при характерных для них значениях плотности 


р =1500-1800 кг/м?, модуля Е =10°-10? Па и скорости ъ, =100-500 м/с, максимально допустимая 
скорость вертикальных деформаций у „, =1-2 м/с [6]. При этом при диаметре излучателя сейс- 
моисточника один метр время действия деформаций лежит в диапазоне от {„, =(1-5):10° с до 


Е ах =(4-20)-10`3 с (середина диапазона (4-7):10`3 с), сила воздействия на грунт ограничена зна- 


чением (2-5)'10° Н, а максимальное значение упругих деформаций грунта 
х к =(2-4)`103 м. 

Следовательно, лишь при малой скорости движения излучателя и малом давлении на по- 
верхность грунта механическая энергия сейсмоисточника независимо от его размеров и мощности 
может использоваться эффективно, так как интенсивность генерируемых волн прямо пропорцио- 
нальна величине упругих деформаций, вызвавших колебания. При более высокой скорости воз- 
действия на грунт в зоне контакта излучателя с ним происходят пластические деформации грунта 
(сжатие под излучателем и сдвиг по его периметру) до глубины, на которой скорость смещения 
частиц грунта уменьшается до величины у „. «Избыточная» механическая энергия сейсмоисточ- 


ника не участвует в создании целевых сейсмических волн, а лишь увеличивает интенсивность 
поверхностных волн — помех — приводит к увеличению запаса «реактивной» механической энер- 
гии в конструкции сейсмоисточника. 

Максимальная длительность воздействия сейсмоисточника ограничена также тем, что его 
излучатель для обеспечения необходимой жесткости при заданных площади контакта с грунтом и 
площади, на которую воздействует импульс силы механического привода, может иметь значи- 
тельную массу. Колебательная система — массивный излучатель на упругом грунте («очаг излуче- 
ния» [4, 5]) — имеет частоту собственных колебаний: 


во = (и, + т,)е, , (5) 
где т, — масса излучателя, кг; т, — «присоединенная масса» грунта [8], кг; е, — податливость 


грунта под излучателем, м/Н. 

Если привод сейсмоисточника создает силу воздействия на излучатель, направленную 
вниз, в течение времени, большего половины периода собственных колебаний этой колебатель- 
ной системы, то часть его механической энергии расходуется не на создание сейсмической волны, 
а на торможение движущегося вверх излучателя, т.е. неэффективно. Величина { соответст- 


тах / 
вующая по формуле (5) длительности одного полупериода: 1„ =л/®, ‚ определяется также по 
формуле (2). Если 1 


облегчать. 


определенная по формуле (2), меньше л/®, , то излучатель необходимо 


тах /! 
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Очаг излучения играет также роль фильтра, подавляющего колебания, частота которых 
отличается от частоты ‹, [4]. Поэтому изменение размеров и формы излучателя в не меньшей 
мере обусловлено техническими требованиями к частотным характеристикам излучаемого сейс- 
моисточником сигнала. 

Известные модели грунта. Для нужд строительства в науке о грунтах используются модели с 
сосредоточенными параметрами, разработанные для определения его несущей способности, в 
том числе и при воздействии на него импульсных или вибрационных нагрузок. Но известные ме- 
тодики определения параметров этих моделей в статике для решения поставленной в исследова- 
нии задачи неприменимы. Можно воспользоваться лишь элементами структуры этих моделей. Все 
такие модели грунта представляют его как систему, обладающую упругими и диссипативными 
свойствами [9]. Некоторые модели просты: модели Кельвина — Фойгта и Максвелла (рис.1, а, 6), 


некоторые сложнее. Например, модель Шведова (рис.1, в), которая содержит кроме элемента е,, 


моделирующего упругие свойства грунта, и элемента 5,, моделирующего вязкое трение в грунте, 


элемент бр, моделирующий пластические изменения грунта. 


5у 


а) 6) | в) 


Рис.1. Статические модели грунта 


Для моделирования процесса взаимодействия излучателя вибрационного сейсмоисточника 
с грунтом используется т„е,5, модель нагрузки (рис.2) [5], где т, — масса прилегающего к из- 
лучателю грунта, колеблющегося в фазе с излучателем, кг; е, — податливость грунта, м/Н; 5, — 
«сопротивление излучения», Н'с/м: 
т, =А„рК” =1,3ЗрЁ*, (6) 
| 12 


е = я, 7 
я рК ъ, р ъ, 7 





5, = К,рЕ”У, = РКУ, ; (8) 


где А -— радиус излучателя, м; А 


т! 


К, , К, — эмпирические коэффициенты. 
Ее Г,С,К, аналог в системе электромеханических аналогий «сила К <> напряжение И », 


представлен в докладе академии наук [10]. 
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Рис.2. Линейная модель грунта под сейсмоисточником 
и ее электрический аналог 


Линейная т,е,5, модель упругого полупространства применима в исследованиях сейс- 


моисточников лишь при соблюдении ограничений (1)-(4). Для более широкого диапазона измене- 
ния параметров воздействий на грунт определение параметров модели грунта (6)-(8) вызывает 
трудности. Продольная волна деформаций, инициированная вертикальным движением излучате- 
ля импульсного сейсмоисточника, при скорости смещения частиц грунта больше у „ [см. форму- 
лу (1)] является в значительной мере волной пластических деформаций. Из-за них возрастает 
масса прилегающего к излучателю грунта т, и может снижаться частота ®, [см. формулу (5)]. 
Лишь после уменьшения скорости смещения частиц грунта на некотором расстоянии от излучате- 


ля до значения у „ продольная волна становится волной в основном упругих деформаций. 


х. 


Предлагаемая модель грунта. Для учета отмеченных свойств колебательной системы «излу- 
чатель сейсмоисточника - грунт» ее расчетную электрическую схему замещения можно составить 
из следующих блоков (рис.3). 





Рис.3. Структура предлагаемой модели грунта 


Блок [СВ выполняет функцию упругой опоры для излучателя сейсмоисточника - ИС. Это 
модель типа т‚е,5, (см. рис.2), формирующая «график колебаний». Присоединенную массу т, 


грунта зададим как нелинейный параметр. Пусть ее величина изменяется от нуля до рассчитанно- 
го по формуле (5) значения за время прохождения волной расстояния, равного диагонали излу- 
чателя [8]: 


2® 2 
при #<—: т, =0,65ру К. 


у 
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Блок $\М/! моделирует с помощью управляемых электрических ключей контакт излучателя 
с грунтом. Излучатель при движении вверх отрывается от грунта в момент, когда сила их взаимо- 
действия становится нулевой, и вновь соприкасается с ним, вернувшись в исходное относительно 
него положение. Для сейсмоисточника, привод ПС которого реализует ударный способ воздейст- 
вия, добавляется показанный пунктиром блок $\\/>, моделирующий контакт и процесс частично- 
упругого соударения падающего груза с излучателем. 

Блок \ моделирует снижение напряжений в сжимаемом грунте до уровня, задаваемого 
пределом упругости грунта с’,. Если скорость ъ, движения излучателя превышает рассчитанную 


по формуле (1) величину, то избыточная скорость «гасится» с помощью функционального источ- 
ника ./ (тока <> скорости [10]). Силу тока такого источника определим функцией вида: 


Л=И (у, >у 


Выходной блок представляет собой линию задержки Т, нагруженную волновым сопротив- 


У, Ух › 0) ) 


тах ? 


лением К, грунта, равным ру, . Сигнал на линию задержки подается с резистора К. блока 1СВ, 


т.е. является пропорциональным скорости смещений частиц грунта. Регулированием затухания в 
линии задержки и ее длины (времени прохождения волны) можно настроить модель на известный 
по результатам полевых сейсморазведочных работ на эталонном профиле случай и в последую- 
щем использовать его в качестве базы для сравнения эффективности одиночных сейсмоисточни- 
ков или их групповой работы. На практике фиксация сейсмоприемниками отражений сейсмиче- 
ских волн с очень большой глубины недоступна из-за наличия сейсмических шумов. Следователь- 
но, по результатам моделирования можно будет считать более эффективным тот сейсмоисточник, 
сигнал с выходного блока которого будет превышать заданный уровень шума по прошествии 
большего времени. 

Моделирование групповой работы сейсмоисточников в этой модели достигается возмож- 
ностью подключения к общему блоку К, выходов линий задержки Т от схем замещения несколь- 
ких сейсмоисточников. 

Заключение. Предложенная модель может быть реализована как в программах численного рас- 
чета систем дифференциальных уравнений, так и с применением программных комплексов элек- 
трических схем, например, Мсго-Сар, Р5рсе и т.д. [11]. В настоящее время продолжается сопос- 
тавление результатов, получаемых с ее использованием, с экспериментальными. 
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